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13  ABSTRACT 

The  nonlinear  development  and  saturation  of  a  single  Langmuir  wave  driven  unstable 
by  a  gentle  bump  in  the  tail  of  the  distribution  function  in  a  collisionless  plasma 
is  studied  by  treating  the  resonant  particles  numerically.  Over  a  wide  range  of  para¬ 
meter  values,  the  airplitude  of  the  potential  $  is  found  to  saturate  at  such  a  level 
that  the  ratio  g  =  u^/yq  =  3.2,  where  w.  «  (ek2$/m)12  is  the  bounce  frequency  of  the 
trapped  particles  in  the  wave  trough  and  yq  *s  the  linear  growth  rate,  approximately 
given  by  the  classical  Landau  value.  In  view  of  the  importance  of  inverse  Landau 
damping  for  many  instabilities,  this  work  should  have  wide  applicability  and  the 
results  should  be  suitable  for  direct  experimental  tests. 
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We  report  here  the  result  of  a  study  of  the  evolution  and  saturation  of 
a  Langmuir  wave  driven  unstable  by  the  resonant  particles  whose  velocities 
are  near  the  wave  phase  velocity  with  a  "gentle  bump-in-tail"  distribution. 
The  method  used  represents  a  new  approach  to  computer  stimulation  of  a 
plasma  in  that  only  a  small  fraction  of  the  particles,  the  resonant  particles 
(i.e.,  trapped  and  nearly  trapped  particles  in  the  nonlinear  stage),  are 
followed  numerically,  the  remainder  being  treated  analytically.  (1)  In  the 
linear  theory,  the  interaction  of  the  wave  with  the  resonant  particles  in 
the  bump  leads  to  the  well-known  Landau  growth  of  the  wave. ^ 

As  the  unstable  wave  grows,  the  nonlinear  effects  become  important.  If 
the  spectrum  of  unstable  waves  is  sufficiently  wide,  then  the  nonlinear 
development  can  be  adequately  described  by  the  quasi-1 inear  theory, (3)  in 
which  the  trapping  of  the  particles  by  the  wave  field  is  neglected.  If,  on 
the  other  hand,  the  spectrum  of  the  unstable  modes  is  narrow,  either  because 
the  linear  growth  rate  is  a  peaked  function  of  wave  number  or  due  to  initial 
conditions,  then  a  single  mode  dominates  and  the  trapping  of  the  resonant 
particles  is  frequently  the  most  important  nonlinear  effect. 

The  effects  of  the  trapped  particles  on  the  damping  of  a  single  large 
amplitude  Langmuir  wave  have  been  previously  studied  with  anal;  tc  methods.  (4) 
In  these  analytic  treatments,  however,  it  is  assumed  that  both  the  amplitude 
and  the  phase  of  the  wave  are  constant  in  time  (yj^  «  1).  With  these 
approximations,  many  aspects  of  the  nonlinear  wave-particle  interaction 
which  are  important  in  the  case  of  unstable  waves,  such  as  the  effect  of 
transition  between  the  trapped  states  and  untrapped  states  (trapping  and 
detrapping)  and  the  adiabatic  heating  and  cooling  of  the  trapped  particles 
due  to  the  amplitude  variation,  are  neglected.  In  the  present  work,  we 
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follow  the  evolution  of  an  unstable  wave  by  treating  the  dynamics  of  the 
resonant  particles  numerically;  it  Is  thus  possible  to  allow  both  the  ampli¬ 
tude  and  the  phase  of  the  wave  to  vary  in  time,  as  dictated  by  Poisson's 
equation,  l.e.  ,  to  study  the  complementary  case  Yg/^g  ^  1  ,  where  is 
ui^  at  t  •  0,  The  nonlinear  phase  shift,  varying  on  the  same  time  scale  as 
the  amplitude,  is  particularly  significant  in  affecting  the  nonlinear 
evolution  and  saturation  as  well  as  in  conserving  the  energy  in  the  wave 
frame.  Neglecting  the  phase  shift  in  the  nonlinear  stage,  as  assumed  in  a 
recent  work/^  leads  to  an  erroneous  conclusion  on  saturation. 

The  bulk  of  the  electrons  (nonresonant  particles)  are  adequately  de¬ 
scribed  by  the  usual  linear  approximation,  even  at  the  time  of  saturation, 
so  they  appear  simply  via  a  dielectric  constant  which  modifies  the  field 
produced  by  the  resonant  particles.  Their  thermal  velocity  is  assumed  to  be 
so  small  compared  with  the  phase  velocity  of  the  wave,  U  *  uWk,  that  they 
can  be  approximated  by  a  cold  fluid  with  dielectric  function 

e  (to, k)  -  1  -  (u )p/oj) 2  .  (1) 


For  the  resonant  particles,  we  assume  a  distribution  uniform  in  space  and 
linear  in  velocity  over  a  range  -6v  <  v-U  =  w  <  5v  : 

(  (a/26v)  (l+6w/5v)  for  |w|  £  6v 


fR(v)  - 


(2) 


0  for  I w |  >  6v 

where  w  is  the  velocity  in  the  initial  wave  frame;  a  is  the  fraction  of 
the  resonant  particles  compared  to  the  bulk  electrons;  8  is  the  fractional 
change  in  f  from  v-U  to  v  •  l  !  jv  ;  and  6v  is  typically  chos«.n  to 
be  of  order  10  (YQ/k)  so  that  there  remain  untrapped  resonant  particles  tn 
the  nonlinear  stage.  This  "single-wave  approximation"  is  valid  for  such  a 
broad  gentle  bump,  where  one  would  expect  many  modes  to  grow,  if  a  single 
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wave  with  amplitude  sufficiently  above  the  thermal  noise  level  Is  launched 

into  the  system  initially,  with  electric  field  K  =  E  sin  kx  and  a  corre- 

o 

sponding  perturbation  in  the  charge  density  of  the  nonresonant  particles: 

Pmd(x)  “  P  (0)  cos  kx  where  p  (0)  =  kb  /4ti  .  The  evolution  of  the  wave  in 
NK  o  o  o 

time  is  described  by  Poisson's  equation 


e(w,k)  E(w,k)  -  47i[po(u),k)  +  pR(ui,k)]/ik  (3) 

where  p  is  the  charge  density  of  the  resonant  particles.  Transforming 
K 

back  to  the  time  domain  and  using  Eq.  (1),  we  have 

-itu  t  -iiu  t  ft  iui  t ' 

E(k,t)  -  ( E  / 2 1 )  e  P  +  (ATT/ik)  [p_(t>  +  (u>  / 2 i)  e  p  dt  e  P  p„  ( t ' )  ] 

O  K  p  K 

o 

(4) 

where  the  solution  with  phase  velocity  -U  has  been  dropped  because  it 
cannot  have  significant  interaction  with  the  resonant  particles  which  have 
v  =  U  .  The  charge  density  of  the  resonant  particles  is  given  by 

N  -ikx  (t) 

PR(t)  =  (anq/N)  J  e  •*  (5) 

j*l 


where  n  is  the  total  particle  density  and  x^(t)  is  the  instantaneous 

position  of  the  resonant  particle.  It  is  convenient  to  express  x^(t) 

in  terms  of  the  deviations,  Zy  from  the  unperturbed  orbit  of  the  j t 

particle  with  initial  velocity  U  +  w,  :  kx ,  *  to  t  +  ill.  +  z,  where  1 1),  = 

J  j  P  J  j  J 

kfx^o  +  w^t)  is  the  unperturbed  phase  in  the  initial  wave  frame.  Since  the 
time  scale  for  the  energy  transfer  between  the  wave  and  resonant  particles, 


typitally  of  the  order  of  the  inverse  growth  rate  or  bounce  frequency,  is 

much  longer  than  the  plasma  period  '0  ^  ,  it  is  convenient  to  remove  the 

P 

explicit  factor  exp(-iw  t)  from  Eqs.  (4),  (5),  by  setting  PD(t)  = 

P  K 

exp(-iiu^t)  p(t)  and  E(k,t)  ■  F(Eo/2I)  exp(-iiu^t)  .  Note  that  F  is  in 
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general  complex,  F(t)  -  |f|  exp[i(j>(t)l  with  F(0)  =  1  .  Choosing  the 
Initial  bounce  time  ,  as  a  unit  of  time,  k*1  as  a  unit  of  distance, 

we  have  the  following  dimensionless  Poisson  equation: 


iuj  t 

(2iE(k,t)/Ej  e  P  5  F(t)  =  1  +  (2ctu/N) 


with 


S(t) 


N  -i<l|)  +z  ) 
I  e  J  i 

j'l 


The  equation  of  motion  is 


[ S(t)  +  (U/2i)  j  dt’  S(t’)] 
0 

(6) 

(7) 


H'p.+z  ) 

Zj  =  (F/2i)  e  J  J  +  c.c. 


(8; 


The  initial  conditions  for  (8)  are  Zj(0)  -  ij(0)  «  0.  To  simulate  the 
distribution  function  of  the  resonant  particles,  Eq.  (2),  we  choose  the  Ip 
to  be  Nz  uniformly  spaced  points  on  the  interval  (0,2tt):  ip^  -  2ni/N  , 

0  <  J,  _<  (Nz~  1)  .  For  each  Ip^  ,  we  choose  values  of  w^  on  the  interval 
(-6v,<5v),  selected  to  correspond  to  the  resonant  particle  distribution 
function,  Eq.  (2): 

_  ...  2  1/9 

(9) 


Wj/Sv  =  [(1  -  1/8)  2  +  4j/6Nv]1/l2  -  1/g 


It  will  be  shown  in  the  following  that  this  distribution  of  N  »  N  N 

z  v 

discrete  particles  does  give  approximately  the  same  linear  growth  rate  and 
frequency  shift  as  the  continuum  distribution,  Eq.  (2). 


Linear  Theory 

For  sufficiently  early  times,  we  can  linearize  not  only  the  nonresonant 
particle  motion  but  also  that  of  the  resonant  particles,  whereupon  (6)  and 
(8)  lead  to  the  usual  multi-beam  dispersion  relation 
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ex(v) 


i  -  (aU3/2Nyv)  ^’(v-Wjf 


(10) 


where  the  sum  is  carried  out  only  over  the  different  values  of  w  given 
by  Eq.  (9). 

In  the  continuum  limit  (Nv  -*■  ») ,  corresponding  to  the  choice  (2)  for 
fR,  the  integration  over  w  in  (10)  gives 


e2(v) 


1  - 


2aU 

v6v 


v8  +  6v 


+  h  log 


v-6v 

v+6v 


(11) 


.  v  -  (6v) ' 

If  6v  «  v  ,  we  expand  the  terms  in  brackets  in  powers  of  <$v/v  , 
obtaining  to  leading  order  the  usual  single  beam  dispersion  relation.  This 
case  has  recently  been  studied  by  O’Neil  et  al. (1)  with  an  approach  similar 
to  ours.  In  the  limit  6v  >>  v  ,  we  can  expand  in  powers  of  v/6v  , 
obtaining 


£,(v)  =  1 


2aU 


V(6v) 


2  {  1  ■  1  I®  +  tr  +  t  <v/6v)3]  } 


The  only  root  of  *  0  is 


v  =  YL(i-2/ir6)  [1  +  4yl/tt6v] 


-1 


(12) 


(13) 


where 


Yl  m\  U2  fQ’(U)  wp  «  toBU3/4(6v)2  (14) 

is  just  the  Landau  value  (in  our  dimensionless  units)  with  f ^ ' (U)  computed 
from  Eq.  (2). 

We  see  from  Eq.  (13)  that  the  growth  rate  is  approximately  the  Landau 
value  provided  «  1  .  The  deviation  of  Yq  from  the  Landau  value  yl 

and  the  nonvanishing  real  part  of  v  are  due  to  the  sharp-edged  distribution, 
Eq.  (2).  This  sharp  edge  can  be  removed  by  joining  smoothly  the  distribu¬ 
tion  of  the  resonant  particles  to  that  of  nonresonant  particles.  This 
introduces  additional  "boundary  terms"  in  Eqs.  (3)  and  (6)  and,  in  the 
linearized  case,  yields  the  exact  Landau  growth  rate  (14),  with  vanishing  Re  v 
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IIh:  near  vanishing  of  at  early  times  In  Fig.  2  1h  a  consequence  of  the 
fact  that  our  numerical  results  are  based  upon  a  'odified  form  of  (6)  which 
includes  these  "boundary  terms".  f  they  are  omitted,  then  a  non-vanishing 
$  ,  of  order  y  ,  is  obtained  in  the  initial  stages,  in  agreement  with  (13). 

Aside  from  this  non-zero  $  and  the  small  changes  in  the  linear  growth  rate, 
however,  those  boundary  terms  do  not  affect  the  nonlinear  saturation  level. 

It  is  determined  by  the  actual  linear  growth  rate  alone,  given  by  either 
Eq.  (13)  or  Eq.  (14),  corresponding  to  the  cases  with  or  without  the  sharp 
edge  in  the  distribution. 

We  return  now  to  the  dispersion  Eq.  (12)  for  "beams".  In  general, 
we  expect  to  find  an  unstable  root,  ,  with  real  part  between  each  adjacent 
pair  of  Wj ,  as  illustrated  in  Fig.  1  for  a  particular  choice  of  parameter 
values.  The  associated  residue  values,  [^'(Vj)]"1  are  also  shown.  The 
continuum  value  for  this  case,  obtained  from  Eq.  (13),  is  indicated;  it  falls 
very  close  to  the  fastest  growing  root  of  the  Nv-beam  modes.  The  latter  also 
has  the  largest  residue  [ e ^ J  \  and  hence  dominates  the  linear  behavior, 
leading  to  a  good  agreement  with  the  continuum  results. 

Nonlinear  Results 

We  put  (6)  and  (8),  including  "boundary  terms",  into  standard  difference 

equation  form  and  solve  by  straightforward  step-ahead  in  time.  The  wave 

initially  grows  at  the  growth  rate  y  given  by  Eq.  (16)  until  a  time  t 

°  NL 

at  which  the  orbit  of  the  particles  trapped  at  the  trough  of  the  wave  becomes 
nonlinear.  The  time  for  the  onset  of  nonlinearity,  t^ ,  is  approximately  given 
j-CNL 

by  J  dt  u>b(t)  =»  or  <ytNL)  ~  71 V 2  +  1  '  After  CNL  ls  reached- 
0 

the  growth  rate  begins  to  decrease  from  the  linear  value,  turns  negative,  and  then 
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osc  illates  about  zero  with  the  instantaneous  bounce  frequency,  as  shown  In 
fig.  2.  The  amplitude  of  the  electric  field  is  initially  small,  i.e., 
g(t-O)  -  UJbc)/Y0  <  1  •  It  Increases  exponentially,  with  growth  rate  v 

o 

untU  tHL,  anc  then  saturates  at  such  a  level  that  the  bounce  frequency  of 
the  trapped  particle  in  the  wave  trough  is  approximately  equal  to  3.2  times 
the  linear  growth  rate,  i.e.,  g  -  ^  .  3.2.  Thls  ratlo  g  at  saturatlon 

is  independent  of  yo  and  other  parametric  values  so  long  as  g  <  1 ,  as 
shown  in  Table  I.  The  nonlinear  frequency  shift  is  of  the  same  order  of 
magnitude  as  the  nonlinear  growth  rate. 

In  conclusion,  we  note  that 

1)  Treating  numerically  only  the  resonant  particles  is  an  efficient 
simulation  technique  for  any  problem  where  a  small  fraction  of  particles 
exchange  energy  with  waves. 

2)  The  universal  character  (g  =  3.2)  of  the  saturation  level  should 
be  suitable  for  experimental  test. 

3)  Since  many  other  instabilities  involve  the  same  wave-particle 
resonance  mechanism  as  Landau's  problem,  similar  results  are  to  be  expected 


there  as  well. 
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Fig.  l 


Fig.  2. 


Figure  Captions 

a)  Locus,  in  the  complex  m  plane,  of  roots,  v^,  of  the  multi- 
beam  dispersion  Eq.  (10)  for  -  30,  with  other  parameters  equal 
to  those  in  Table  1,  Case  A.  The  vertical  lines  show  the  locations 
of  the  30  beams.  The  point  designated  with  an  arrow  is  the 
continuum  value  of  v  given  by  (13). 

b)  Relative  values  of  the  residues,  for  the  roots 

of  (1)  shown  In  a).  As  in  a),  the  second  highest  point  is  the 
residue  obtained  from  the  continuum  form,  (12),  for  c. 

Temporal  evolution  of  the  instantaneous  frequency,  uj  (t)  « 

b 

1/2 

[eE(t)k/m]  ;  growth  rate  y  ,  and  frequency  shift,  $  ,  for 
case  A  of  Table  I.  with  -  960  and  Nz  -  4.  The  vertical  scale 
is  in  the  unit  of  initial  bounce  frequency,  iu.  ;  the  horizontal 

DO 

scale  is  in  the  unit  of  .  The  Landau  value,  y^,  given  by 

(14)  is  also  shown. 


Ac  Summary  of  parameters  and  saturation  amplitudes  for  various  cases 
1/2 

^bo  "  (ekhQ/m)  is  the  initial  bounce  frequency,  is  the 
bounce  frequency  at  saturation,  yL  is  the  Landau  growth  rate, 

Yq  the  actual  initial  growth  rate. 
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Fig.  1.  a)  Locus,  in  the  complex  w  plane,  of  roots,  v^,  of  the  multi-beam 

dispersion  Eq.  (10)  for  =  30,  with  other  parameters  equal  to  those  in  Table  I, 
Case  A.  The  vertical  lines  show  the  locations  of  the  30  beams.  The  circled 
point  designated  with  an  arrow,  is  the  continuum  value  of  v  given  by  (13). 
b)  Relative  values  of  the  residues,  [e  ^ ’  (v  )  J  1  for  the  roots  of  (1)  shown  in 
a).  As  in  a),  the  second  highest  point  is  the  residue  obtained  from  the  con¬ 
tinuum  form,  (12),  for  c. 


